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In this paper we develop a mathematical framework for 
the characterisation of separability and entanglement of for- 
mation (EoF) of general bipartite states. These characteri- 
sations are of the variational kind, meaning that separability 
and EoF are given in terms of a function which is to be mini- 
mized over the manifold of unitary matrices. A major benefit 
of such a characterisation is that it directly leads to a numer- 
ical procedure for calculating EoF. We present an efficient 
minimisation algorithm and an apply it to the bound entan- 
gled 3x3 Horodecki states; we show that their EoF is very low 
and that their distance to the set of separable states is also 
very low. Within the same variational framework we rephrase 
the results by Wootters (W. Wootters, Phys. Rev. Lett. 80, 
2245 (1998)) on EoF for 2x2 states and present progress in 
generalising these results to higher dimensional systems. 

03.65.Bz, 03.67.-a, 89.70.-l-c 



I. INTRODUCTION 

A problem which has received considerable attention 
in the last few years is to find necessary and sufficient 
conditions under which a quantum state of a compos- 
ite system is separable. The example extraordinaire of a 
non-separable state is a pair of 2-level particles in a sin- 
glet state, a so-called EPR-pair, named after Einstein, 
Podolsky and Rosen, who used this sort of state to show 
that quantum mechanics exhibits strong non-local corre- 
lations which seem to violate the relativity principle. 

A separable state of a composite system can be written 
as the direct product of the subsystem states: I^'ab) = 
l^*^) ®\^b)- a non-separable state, or entangled state 
cannot be decomposed in this way; e.g., the singlet state 
(I T)l i) ^ I i)l T))/%/2 consists of a superposition of 
separable states but is itself not separable. 

Nowadays, the importance of entangled states goes be- 
yond a mere fundamental interest, since EPR-pairs are 
the basic resources of quantum techniques such as quan- 
tum cryptography, quantum teleportation and quantum 
error correction. A mixed state is separable iff its density 
matrix can be written as a convex linear combination of 
pure product states; for a bipartite system this reads: 



K 
k=l 



Wk\u''){u''\(g>\v^){v''\. 



(1) 



The separability problem consists of finding a criterion 
for checking whether such a decomposition is possible for 
a given state. 

Despite the simple formulation of this problem, a com- 
plete solution has to this date not been found. An im- 
portant achievement was the discovery by Peres of a nec- 
essary condition for separability He noted that the 
partial transposition of a separable state still has non- 
negative eigenvalues, just as the original state. Thus, if 
the partial transposition of state p is not a state (i.e. does 
not have non-negative eigenvalues summing to one) ,then 
p is not separable (i.e. is an entangled state). The im- 
portance of this criterion was soon realised when M., P. 
and R. Horodecki proved ||^ that it is also a sufficient 
criterion for 2x2 and 2x3 systems. For an introduction 
to recent results on this subject, see, e.g., 

If a state is entangled, one could ask for a measure of 
the amount of entanglement. For pure states, a measure 
generally agreed upon is the von Neumann subsystem 
entropy: the entropy of the partial trace of the state 
projector. For mixed states, the situation is much more 
difficult. Not only is there no single measure of entangle- 
ment which is suited for every purpose, but calculating 
the values of the different proposed measures and proving 
statements about them is exceedingly difficult. Among 
the proposed measures are the entanglement of forma- 
tion , the entanglement of distillation Q and relative 
entropy of entanglement 

In this paper, we focus on separability, on entangle- 
ment of formation (EoF) and on the related concept 
of concurrence. All these subjects are related, because 
states are separable if and only if their EoF is zero. A 
closed form expression exists for the EoF of 2 x 2 systems 
in terms of their concurrence . A closed form expres- 
sion also exists for isotropic states of general systems . 

The purpose of this paper is to give variational char- 
acterisations of separability and EoF for general (i.e. any 
dimensions) bipartite states. Such a characterisation is of 
the form Q{p) ~ min^ f(p, T), that is: the state property 
under study can be found as the minimal value of some 
specific function over the manifold of unitary matrices T. 
In section II it will be shown how this can be done. The 
language of section II is matrix analysis, not only because 
this allows to state the results in a most succinct way but 
also because it gives clues towards generalisations. 

The greatest benefit of a variational characterisation is 
that it directly yields a method for actually calculating 
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the state property Q, albeit in a numerical fashion, us- 
ing a minimisation procedure. In section III we describe 
the procedure which we have used, and some interesting 
results we have obtained with it. 



circumstances, this method directly yields an optimal re- 
alisation, without any need for searching. We have not 
yet investigated whether this method is applicable to the 
EoF case. 



II. VARIATIONAL CHARACTERISATIONS 



A. Relation between different realisations of a state 



It is well known that mixed states can be realised by 
an ensemble of pure states in an infinite number of ways. 
The determination of the separability of a state and the 
determination of its entanglement of formation have in 
common that a particular realisation of a state has to be 
found such that some property holds for all pure states in 
that realisation. In order to find this optimal realisation, 
it is of considerable interest to have a mathematically el- 
egant way of "g enera ting" all possible realisations of a 
state. In section [I A we will recollect a result by Hugh- 
ston, Jozsa and Wootters that any realisation of a state 
is related to the eigenvalue decomposition of the state via 
some right-unitary matrix. 

The required property for separability is that all pure 
stat es in the realisation must be product states. In sec- 
tion [I B we give a number of useful mathematical expres- 



sions for this property. This then leads to a variat iona l 
characterisation of separability, the topic of section [I C . 
For calculating the EoF of the state, the property of the 
optimal realisation is that the so-called average entangle- 
ment of the realisation is minimal. This property and an 
ensuing variational characterisation of EoF will be dis- 
cussed in section 
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In this way, searching all possible realisations for some 
property amounts to passing through all right-unitary 
matrices and test the property in question. However, 
this would be a very impractical way to determine sepa- 
rability or EoF if there would not be some bound on the 
dimension of these matrices, or, which is the same thing, 
on the number of pure states in the optimal realisation. 
Luckily, such a bound exists. In the case of separabil- 
ity, Horodecki proved that (7ViiV2)^ pure states (or 
less) suffice, where Ni and A^2 are the dimensions of the 
subsystem Hilbert spaces. Uhlmann proved that a 
similar bound holds for the determination of EoF: the 
number of pure states in the optimal realisation need not 
be larger than the square of the rank of the state. 

In section II E we discuss the so-called concurrence of a 



state, a quantity which is closely related to the EoF. We 
give an alternative proof of an important result on the 
concurrence of 2 x 2 states by Wootters Q . One of the 
virtues of this alternative proof is that it yields an addi- 
tional result on the exact amount of pure states in the 
optimal realisation. We then report some progress in gen- 
eralising the concurrence concept to higher-dimensional 
bipartite states. 

In appendix A, finally, a method is described for re- 
ducing the set of unitary matrices which has to be ex- 
amined, in the case of separability testing. Under some 



Consider a rank-i? state p in an A^i x iV2-dimensional 
Hilbert space, realised by an ensemble {wk,\tp'^)})^^i, 
where the Wk are the mixing weights of the K pure state 
vectors I'tp'^). The number K is called the cardinality of 
the ensemble. Necessarily, K cannot be smaller than the 
rank R. Since there generally are an infinite number of 
ensembles realising a particular mixed state, we are free 
to choose K larger than R if this suits our purposes. It 
will turn out that sometimes we will even be forced to 
take K > R. 

Thus: p = Y.k=iWk\ij''){ip''\, or p = ^W^'' , where 
W is a. K X K diagonal matrix with Wkk = Wk, and the 
columns of ^' are the K vectors ip''. This decomposi- 
tion of p is reminiscent of the eigenvalue decomposition 
of p: p = ^M^t^ where M is an i? x i? diagonal ma- 
trix whose diagonal elements are the eigenvalues of p, 
and the columns of are the R eigenvectors. Since p is 
Hermitian, $ is a unitary matrix. 

It can now easily be proven that these two decomposi- 
tions must be related by an i? x if right- unitary matrix T; 
this has been done first by Hughston, Jozsa and Wootters 
§. 

Lemma 1 For a general state p, with eigenvalue decom- 
position p — $A/<I>t, there is a matrix ^ and a non- 
negative diagonal matrix W such that p = ^'VF^'^ iff 
there is an Rx K matrix T such that: 



^1^1/2 ^ $M^/^T, with TT^ = Ij 



(2) 



Right-unitarity of the matrix T means that a unitary 
K X K matrix T' exists such that T consists of R row 
vectors of T; that is, the R row vectors of T form an 
orthonormal set in and the K column vectors are 
projections of an orthonormal basis in onto an i?- 
dimensional subspace. Stated in matrix algebraic terms, 
the proof becomes very simple: 

Proof. First of all, it is obvious that <i>Af<i>1' = vj/jyvl/t 
follows directly from (|^). Conversely, denote X — vj/iy^/^ 
and consider the singular value decomposition oi X: X — 
C/Sy, where J7 is a unitary R x R matrix, V a right- 
unitary R X K matrix and E a diagonal R x R matrix 
with non-negative diagonal elements. From ^M^^ = 
^-Vt^^ft we get $M$t UY?W . Since both M and S are 
positive semidefinite, it follows that S = U^^M^/'^'^^U 
so that X = ^M'^/'^^^UV. This is precisely equation(||), 
with T = □ 

Remark. It is noteworthy that the elements of W and 
M are related to each other independently of and 4": 

Wk = {T^MT)kk- 
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This follows from the observations that <I> is unitary and 
that the columns of '5 have norm one. 



B. Characterisation of product states 

A state of an iVi x N2 system is separable iff there 
exists a realising ensemble consisting solely of product 
vectors ^ — with tp-^ g Hi and -0^ S "^2 (in 

this paper we use superscripts for enumerating vectors, 
and subscripts for denoting vector components). Product 
vectors can be characterised easily by rearranging their 
components in matrix form. For an 7ViA^2-vector x, let x 



be an A^i x 7V2 matrix such that x 
product vectors this gives: 



Yl^ij iy-e* <Si e^. For 



Obviously, product vectors are characterised by the con- 
dition that the rank of ■(/' is 1. A necessary and suf- 
ficient condition for this is that all 2 x 2 minors of ij; 
must be zero, or, more succinctly, that the second com- 
pound matrix of i/j must be zero: C2('0) = 0. The 
second compound matrix of an A^i x iV2 matrix is an 
{Ni{Ni ~ l)/2) X (iV2(iV2 - l)/2) matrix with elements: 



C2(^)(. 



A-tj'Ai'j, 



The elements of C2 are all possible 2x2 minors of 
A. The second compound matrix has a lot of useful 
properties, such as: C2{AB) = C2{A)C2{B), C2(l„) = 
K(n-i)/2 and C2{A-^) = {C2{A))-^ §. 

For practical applications it is sometimes better to con- 
sider a {Ni — 1) X {N2 — 1) submatrix of C2, the one con- 
taining the elements C2(A)(i q j+i) only. It is easily 
seen that the vanishing of this submatrix is already suf- 
ficient for A being of rank 1. 

From the expression for the second compound matrix, 
which is quadratic in ^, it will prove useful to construct 
a bilinear function of two A^i x N2 matrices, denoted 
C{A,By. 



C{A,B)(^,i').{j]') 



Aij,Bi,j, i <i',j < f 



Obviously, C{A,A) = C2(A), so that C{A,A) = if and 
only if A has rank 1. More specifically, we can apply 
this to the state vectors tp'^: tp'^ is a product vector iff 
Cii}'',^!'') = 0. 

In the following, we will only use a symmetrised version 
of C, which we will denote by 

Since this is a bilinear function in the elements of we 
can express this in matrix notation: 

where the notation (a) is a shorthand for the index tuple 
The matrices S*^"-*, which we call indicator 
matrices, are defined as 



c.(") 



= 5, 

= S. 



(a) 



(a) 

(i'j)Xii') 



= -1 



all other elements being zero. Note that all S have rank 
equal to 4. For the case of 2 x 2-systems, there is only one 
indicator matrix; it is equal to Uy® ay, corresponding to 
a spin-fiip operator [0. 



C. Condition for separability 

We can now formulate a general necessary and suf- 
ficient condition for the separability of a mixed state. 
As mentioned before, the state p = ^M^^ is separa- 
ble iff there exists a decomposition p — '^W'^\ with 
^fVt^i/s = $Mi/2y^ such that aU ip^ are product states, 
or C(?/;'=,7/;') = 0, for aU k^l. 

Now: 

= J2 ^^ci^,r,v^^''), (3) 



p,g 



where we have used bilinearity of the form C. 
Given the eigenvalue decomposition of p, the entity 
C{y/rn^(pP , y^rn^cj)'^) can be calculated straightforwardly. 

Let us organise its components into a set of matrices A*-"' : 



= C{^(pP, V77I^(/)«)„ = ^/M$^5(")$VM. (4) 
Using this notation, (^) can be written concisely as 

C{^\^') = {T^A^'''>T)ik/V^i- 

The state is therefore separable iff we can find an i? x if 
matrix T, with K > R, such that 



(5) 



Here, k ranges from 1 to K, and a enumerates all tuples 
of indices {i,i',j,j') with 1 < i < i' < Ni and 1 < j < 
j' < -^2- As noted before, it is also sufficient to consider 
only the tuples + + 1). 

Testing separability requires that the system (^) be 
solved for T. Another approach, however, is to consider 
(T'^ A^°''^T)kk as entries of a matrix indexed by a and k 
and to try to minimise a matrix norm of this matrix. The 
state is then separable iff this minimum is zero. One can 
use whatever matrix norm one prefers, e.g. the Hilbert- 
Schmidt norm (also called Frobenius norm or ^2-norm) 



\M\l = Y.^i l^ijP = TvAA^ Thus p is separable iff 



mm^|(T^A(")T)fe,p = 0, 

OL.k 



(6) 
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where the minimum has to be taken over all if > i? and 
all Rx K matrices T for which TT^ = Ir. The minimal 
K is called the cardinality of the state. 

One can also use the h norm (sum of absolute values) 
and minimise J2a k \{T'^A^"^T)kk\- For 2x2 systems the 
li norm is the average concurrence of the ensemble, as 
introduced by Wootters in ^J, and the minimum is the 
concurrence of the state p. Note that in the context of 
separability testing it does not matter whether one uses 
(T^^(")T)fcfc or {T^A(-)T)kk/wk. 

To end this paragraph, we derive an alternative expres- 
sion for the h norm ||(C2(?/''^))fc||2- Define B'^ = 'ij)^ 
with eigenvalue decomposition = U^Ti^U^'^ (with 
T,^ — Diag((Tf )). Using the properties of C2 we find 

ll(C2(^'=))Lill^=ETr(C2(^'=)C2(V;'=)t) 

= ^ Tr ^2(5'=) = ^ Tr ^2(2'=) 

k k 

k i,j i 

= ^E(E-f)^-E(-^)^ 

A; i i 

= i^(TrS'^)2-Tr(E'=)2 

k 

= B^f - TY[B^f . 



This result can be interpreted easily: a positive definite 
hermitian matrix is rank 1 iff the square of its trace equals 
the trace of its square. 



D. Entanglement of formation 

Within the same framework, we can also give a varia- 
tional characterisation of the entanglement of formation 
E{p) (EoF) of a mixed state p. This quantity is defined 
as the average entanglement of the pure states in a real- 
ising ensemble, minimised over all possible realising en- 
sembles Q. The von Neumann entropy H of a state p is 
— Trplog2p; introducing the function h{x) = — xlogjX, 
we can express H as a function of the eigenvalues Xk of 
p: H{p) = J2k ^i^k)- The entanglement of a pure state 
■0 of a bipartite system {A, B) is the entropy of the par- 
tial trace of the projector of i?(V') = H{pa), with 
PA = TrB(|'0)(^/'|). The average entanglement of an en- 
semble {wk,^^} is J2k''^kE{ilj''); the EoF is then found 
as the minimal value over all ensembles realising p. 

In this paragraph, we will derive an expression for E{p) 
which is better suited for calculation. Let {wfc,?/'''} be 
the realising ensemble with least average entanglement 



and {rUp, (jjp} the realising ensemble corresponding to the 
eigenvalue decomposition of p. We first express the par- 
tial trace of the projector of ip'' in terms of ip'^: ip'^ = 
Ei,,^'e' « hence = Er.,.,.J'MLne' » 



e^)(e^ (E) e'^y , and the partial trace equals 
p\^TrBm{^'\) 

i,P \ 1 / 

This is precisely the matrix S*"' from the previous para- 
graph. 

Remark: The entropy of this partial trace matrix 
can be expressed in terms of the singular values of tp'^. 
Let Tp'' = U'^'S'^V'^ be the singular value decomposi- 
tion of tp'^ (that is, the Schmidt decomposition of tp'')^ 
with [/'^ unitary and V'^ right-unitary (supposing that 
A^i < N2) and S'' a positive semidefinite diagonal ma- 
trix, then p'x = U'^iY^'^ fiU'^Y and H{p'X) = H{{Y.'' f) = 
-2E,(af)^log2(af). 

In the present framework only the eigenvectors (jf are 
known, and the vectors ip^ are to be sought by looking for 
an appropriate T-matrix. We therefore want to express 
p\ in terms of T and the (pP . We get: 

Wkp\ = y/w^i!^ ^/w^{'|p^)^ 

R 

= E ^Pi'^qk ^/rUpniqipP (0«) . 

p,q=l 

Let us use the symbol Aj. (T) as a shorthand for the right- 
hand side of the previous expression: 



Afc(r) 



R 



p'X = Ak{T)/wk 



Wk = Tr Afc(T). 

The last equation follows from the fact that p\ is nor- 
malised. 

The EoF is thus: 



K 



' k=l 
K 

= mi^EG(Afc(T)), 



fc=i 



with 



G(A) = -Tr(Alog2(A/Tr(A))) 
= H{A) - ft(Tr(A)). 



(7) 



(8) 



4 



The minimum has to be taken over all if > i? and all 
R X K matrices T for which TT^ = 1r. Note that, 
since a state is separable iff its entropy of formation is 
zero, equation (|^) gives an alternative for equation (|^) 
for testing separability. 

Equation (|^) can be brought in a more suitable form if 
we enlarge the set of matrices 0^ with zero matrices for 
p > R. Then we can always use square, and therefore 
unitary T matrices. Following a result by Uhlmann plf , 
the cardinality K must lie between the rank R and the 
square of the rank. This guarantees that the EoF can be 
found by restricting oneself to finite sized T matrices. 

E. Concurrence 



The following theorem by Thompson gives a precise 
relationship between the moduli of the diagonal elements 
of a complex square symmetric matrix and its singular 
values [12| (stated without proof): 

Theorem 1 (Thompson) Let di,...,dn be complex 
numbers and si,...,s„ nonnegative real numbers, enu- 
merated so that \di\ > ■■■ > \dn\ and si > ■■■ > s„. 
A complex symmetric matrix exists with di, . . . ,dn as its 
diagonal elements and si, . . . , s„ as its singular values, if 
and only if 

k k 
i=l i=l 



The first analytic formula for calculating EoF has been 
found by Wootters and is valid for 2x2 systems. A ba- 
sic property used in deriving the formula is the so-called 
concurrence of a state. The concurrence is also useful for 
testing separability, because a 2 x 2 state is separable iff 
its concurrence equals zero. In this section we do two 
things: first we rederive Wootters' results in a shorter 
way, based on the concepts we have introduced above 
and using an interesting theorem from matrix analysis. 
This rederivation gives hints toward the generalisation of 
the concurrence concept to higher-dimensional systems, 
which is the second topic of this section. 



1. r/ie 2 X 2 case 

In this paragraph we give a shorter proof of Wootters' 
results on the EoF of 2 x 2 systems Q. For the case of 
2x2 systems, formula (|]) becomes particularly simple, 
since there is only one 2x2 minor to consider, so that 
there is just a single symmetric matrix A^"^. 

The concurrence of a pure state equals C(V') = 
IV'^'S'V'I- For the pure states ■0'^ in a decom- 



(*^5*) 



fcfc 



position of p, we get C{tp'^) 
\{W-^'^T^ATW^/^)uk\ = \{TTAT)kk\/wu. 

The average concurrence of a realisation of p is thus 
given by '^k\{T'^ AT)kk\ and the concurrence of p is 
the minimal average concurrence over all possible real- 
isations, i.e. over all possible right-unitary T . Since A 
is symmetric, its singular value decomposition assumes a 
special form, known as the Takagi eigenvalue decomposi- 
tion 1^: A = U'^T,U (again, U is unitary and E positive 
semidefinite diagonal). Since we consider all possible T, 
the matrix U can be absorbed in T, so that the expres- 
sion for the concurrence becomes nLinrX^j, \{T'^T,T)kk\- 
So, T^J^T runs through all possible complex symmetric 
K X K matrices with R prescribed singular values S (if 
K > R then K — R zero singular values have to be added 
to S) and the average concurrence equals the sum of the 
moduli of the diagonal elements. 



i—k 



n-3 

E 

i=l 



n~2 



\di \ - \dn^2 \ - \dn~l \ - \dn\ < (^ Si) - S„_i - S„. 



The last inequality does not apply when n < 3. 
The second inequality gives, for k — 1: 



^ Mil > si - (V" Si). 



1=1 



i=2 



Applied to the problem at hand, we find that the minimal 
average concurrence must be fi — (^,^2 O'" zero if this 
quantity is negative. Here we have put K ~ A. Letting 
K be larger than 4 can give no improvement, since this 
amounts to just adding K — A zero singular values, and 
this does not influence the inequalities of the theorem. 

If i? < 4, we could try to put K — but then the 
third inequality comes into play: 

3 3 

i=l i=2 



so that 



Cip) 



K=3 



\<Tl — (72 — CTsI 



Therefore, if i? = 3 and cti — (12 — era < 0, putting K = 
4 gives zero EoF, while K = 3 gives non-zero EoF. In 
other words, these states are separable in (at least) four 
product states {K = 4). Furthermore, a rank 3 state is 
separable in three product states {K = 3) iff cti — (T2 — 
tT3 =0. 

If i? = 2, we can safely put K = 2, since then the third 
inequality does not apply. 
We have thus proven: 
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Theorem 2 The concurrence of a 2 x 2 state equals: 

R 

C{p) = max(0, (Ti - (Ti j). 

1=2 

where Ui are the singular values of the corresponding A- 
matrix, in descending order. The optimal cardinality K 
equals the rank R, except in the case when R — 3 and 
a\ < (72 + era, where the optimal K is 4- 

Because of the statement about the optimal cardinality, 
this theorem is an improvement over Wootter's theorem. 



2. Relation between concurrence and entanglement of 
formation 

For the sake of completeness, we rephrase the rest of 
Wootters' results of in the present setting. 

The entanglement of a pure state is a convex, 
monotonous function £ of the concurrence of the state: 
E{ip) = £{C{ip)). Hence, the EoF, which is the average 
pure state entanglement, equals 



E{p) 



mm > 

k 



Wk£{\iT'^AT)kk\lwu). 



Because of the convexity of f , this gives E{p) > 
mm J' 

£{J2k \iT'^AT)kk\), where equality holds only if all 
quantities KT"^ AT)kk\/wk are equal. Using Thompson's 
theorem again and the monotonicity off, this minimum 
is equal to £{ai - J2j>i '^i) = ^{Cip)). 

We therefore look for an optimal T matrix, yielding 
minimal average concurrence (C(p)), and for which, ad- 
ditionally, all the quantities \{T'^ AT)kk\/wk are equal 
(and thus equal to C(p)). There exists a T' for which 
J2kiT'^^T')kk is equal to C(p); indeed, with A = 
C/^EC/, set UT' = Diag{l, i,i,...,i), then T'^AT' = 
Diag(l, — 1, — 1, . . . , — and the trace of this matrix 
is (Ti — ((72 -t- • • ■ -t- (Jk)- If this quantity is positive, it is 
equal to C(p); if not, p is separable and we immediately 
have that all \{T''^AT')kk\/wk are equal (zero). 

Concerning the non-separable states: for any or- 
thogonal matrix O, Tr{T'0)^ A{T'0) = Tr T'^ AT'. 
As described in using a suitable O we can make 
all {{T'O)^ AT'0)kk equal to a constant a times 
Wk (exploiting the fact that T''^AT' is a real di- 
agonal matrix here). Summing over k then yields 
C{p) = \EkiiT'OrAT'Ohk\ = laEfc^fcl = 
so that ({T'0)^AT'0)kk = C{p)wk. Then, 
\{{T'OYAT'0)kk\ = C(p), so that T = T'O is the 
matrix we were looking for. 



3. Generalised concurrence 

According to equation (^, a state is separable iff a 
right-unitary T can be found such that the diagonal ele- 
ments of every T'^A^^^T are zero. In analogy with defin- 
ing the average concurrence of a realisation of a 2 x 2 state 



as the /i-norm of the diagonal elements of T^AT, in the 
general case we can define a concurrence vector as the 
vector of ^i-norms of the diagonal elements of T^ A^°''^T: 



C(a)(T)=^|(T^A(")T)fc,| 



(9) 



A state is therefore separable iff a T exist such that 
the concurrence vector is zero. From the previous para- 
graph, a necessary condition follows immediately: 



1=2 



V(a), 



(10) 



where the a'f'^ are the singular values of A^°'\ in descend- 
ing order. 

Unfortunately, this condition is not a sufficient one for 
separability. Numerical experiments showed that crite- 
rion (p^ ) is weaker than the Peres criterion, which is a 
non-sufficient criterion itself. The main reason for this 
failure is that all the components of the vector concur- 
rence (^ must be made zero by one and the same T. Typ- 
ically, however, the matrices A*^"^ all have different sin- 
gular vectors (the rows of the U matrix), so that the C/^"-* 
matrices in the decomposition A^") — f/^"^-'"!]*^"''?/^"^ 
cannot all be absorbed in T at the same time. 

It is easy, however, to find a stronger criterion than 
criterion (^0|): as equation (||) is linear in the matrices 
A^"-\ the condition ( p^ ) must hold also for every linear 
combination of the matrices A*-"^ . Denoting the j-th sin- 
gular value (descending order) of the linear combination 
X](q) ^{a)-^''"^ by Uj{x), it follows that another, and po- 
tentially stronger, necessary condition for separability is 
given by: 



max 



0-1 (x) 



a J (x) 



< 1, 



(11) 



where M is the number of tuples (a). Again, one could 
choose to consider all possible A^"^ or just the minimal 
subset with (a) — (i,i + + 1). 

Numerical experiments now showed that criterion ( pl| ) 
is actually stronger than the Peres criterion, provided 
all A*^"'' are used. In the next section we will give an 
example where condition (|l^) even seems to be sufficient 
for determining separability. 



III. NUMERICAL RESULTS 

In this section we present an application of the vari- 
ational characterisations of separability and EoF. Since 
these characterisations involve looking for the minimum 
of a function over a finite-dimensional manifold, it must 
be possible to find a numerical algorithm that actually 
calculates that minimum. As a result, it must be possible 
to calculate the EoF for any bipartite state and, more- 
over, to give the optimal realisation of the state (from 
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the optimal T matrix). In the following paragraphs, we 
first present in some detail a practical minimisation al- 
gorithm for this problem, and then apply the algorithm 
to the calculation of EoF for a family of 3 x 3 states. 



A. Algorithm for minimisation 

Our algorithm for calculating the entanglement of for- 
mation is based on a modified conjugate gradient minimi- 
sation procedure. Starting from an initial point T = Tq, 
conjugate gradient algorithms iteratively seek a direc- 
tion along which progress in minimising the objective 
function is optimal and then perform a so-called line 
search to actually find the minimum along that direc- 
tion. In the present case, however, minimisation is over 
the unitary manifold. This manifold is not Euclidean, 
and the standard line search has to be replaced by a 
geodesic search Q. A geodesic on the unitary mani- 
fold is a one-parameter subgroup of the unitary group: 
T{t) = Toexp(tX), with X a skew-Hermitian matrix 
giving the direction (tangent vector) of the geodesic. 
Through a geodesic search one looks for the optimal t 
for which g{Toexp{tX)) is minimal. 

In steepest descent minimisation, the direction for the 
line search is taken to be minus the gradient of the ob- 
jective function in the current point. Conjugate gradient 
methods improve on this by taking the direction of the 
previous step also in account; if not, the progress made 
in the previous step could be partly undone by the new 
iteration. We have used a modification of the Polak- 
Ribiere formula for calculating the search direction [ p^ ; 
the search direction for iteration i is based on the gradi- 
ent at the current point and on the search direction for 
the previous iteration i — 1: 



X, 
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((Vg).-T(Vg)._i,(V.9).) 



((V5)._i,(V5).-i) ' 

where (,) is the inner product of the embedding space, 
being in this case the standard Hilbert-Schmidt inner 
product (a;,y) — Tra;?/^. The symbol r denotes paral- 
lel transport of the gradient vector at the {i — l)th point 
to the ith point along the geodesic Q : 

T(Vg)i_i = exp(Xj_iii_i/2)(Vg)i_i exp(-Xi_iii_i/2). 

For the line search, we have used the method described 
in again modified to take into account that the 

search is performed along the geodesic g{Tiexp{tXi)). 

Any minimisation algorithm actually finds local min- 
ima. To find the global minimum, we select a number of 
starting points at random and let the minimisation al- 
gorithm work from these points. The minimum is then 
taken over all the results. While this procedure does not 
guarantee that the global minimum is actually found, we 
found that trying about ten starting points gives satis- 
factory results. 



B. Calculation of the gradient 

In this paragraph we give an analytic expression for the 
gradient of the target function g{T). Conjugate gradient 
methods perform better if an explicit expression is given; 
in the absence of such an expression, the gradient has to 
be approximated numerically. 

To calculate the gradient, we have to select an arbi- 
trary direction or tangent vector X, which for the uni- 
tary manifold is a skew-Hermitian matrix. The geodesic 
on the unitary manifold along this direction and passing 
through To is given by T, = Tq exp(eA), or To(l + eX), 
for small e. The gradient of a scalar function on the man- 
ifold can be calculated from the variation of the function 
along the geodesic, using 



de 



where (, ) is the Hilbert-Schmidt inner product. 
The gradient of the target function g{T) is: 



Lemma 2 



where 



and 



g(B,A) = -TrSlog2 



A 
TrA' 



The details of the calculation are given in appendix B. 



C. Results 

As a preliminary test, we have calculated the entan- 
glement of formation of several states of a 2 x 2 system, 
and compared the numerical values with the ones obtain- 
able from Wootters' formula. Furthermore, we consid- 
ered a one-parameter family of 3 x 3 states called isotropic 
states, and compared the numerical values with the EoF 
calculated from Terhal and VoUbrecht's formula 0. In 
all cases, agreement was complete within numerical ma- 
chine precision, except for some isotropic states where 
there was a very small deviation from the formula for 
parameter values close to 8/9. This can be explained by 
the fact that, for these parameter values, there are two 
local minima of the target function which are very close 
in value, and that the minimum with lowest value has a 
very small "basin of attraction" . 

The first interesting results were obtained on the 
Horodecki 3x3 states . These states were introduced 
to show that the Peres criterion is not sufficient for de- 
termining separability. These states exhibit bound en- 
tanglement: their entanglement of formation is non-zero, 
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while their entanglement of distillation is zero (they have 
positive partial transposition). The density matrix of a 
Horodecki 3x3 state is 



p{a) 



1 



1 + 8a 



a 











a 











a 





a 





























a 
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b 



where a is a parameter between and 1, inclusively, and 
b = (1 + a)/2 and c = Vl — a? /2. Note that, since these 
states are not full-rank (their rank is 7), and neither is 
their partial transpose, these states lie on the boundary 
of the set of states and also on the boundary of the set 
of bound entangled states. 

The result of the calculation is shown in Fig. l|. Here 
the entanglement of formation has been calculated of 
a mixture of the Horodecki states with the maximally 
mixed state: ep(a) + (1 — e)l/9. In Fig. |], the scale is 
linear, while in Fig. ^ the scale is logarithmic, so that 
the borderline of the set of separable states is clearly vis- 
ible. The "floor" in the logarithmic picture at -10 is an 
artifact; the algorithm stops when the entanglement gets 
below 10^1". 

Note from these results that the Horodecki states have 
a rather low entanglement of formation (about 0.0109 
for a — 0.225) and that their distance to the manifold 
of separable states is also small (e = 0.93 for a = 0.225; 
that is: mixing the state with just 7% of the identity 
destroys all entanglement). At first sight, the fact that 
the appearance of the set of separable states is not convex 
might seem confusing. However, the parameter a appears 
in a non-linear way in the density matrix so that the 
matrices lie on a non-rectilinear curve in the Euclidean 
state space. The figure, on the other hand, has a as 
parameter and therefore gives a distorted image. 

Fig. ^ shows the entanglement of formation for the 
particular value of a = 0.225 and for e going to 1. From 
this figure, we are led to conjecture that the derivative 
to e becomes infinite at e = 1. 

The abovementioned calculations have been performed 
with the cardinality K set to 14. Fig. ^ shows the effect 
of using different K in the calculations; here e — 1 and 
a = 0.225. It is seen that the value X = 14 is optimal for 
calculating the entanglement of formation in this case. 

For these same Horodecki states, we have also tested 
the conjectured condition for separability (equation 
(|lT|)), based on the generalised concurrence. It turned 
out, quite surprisingly, that the condition correctly pin- 
pointed all separable states, which was verified by com- 
paring the results to Fig. |2[ This leads us to hope that 
equation ( p| ) might be an important step towards find- 
ing a simple and efficient operational criterion for testing 
separability. 



IV. CONCLUSIONS 

We have presented a matrix analytical framework 
within which the questions of separability of mixed states 
and calculating their entanglement of formation can be 
formulated in an elegant and practical way. A main re- 
sult is that, at least in principle, it is now be possible to 
calculate the EoF of any state, or determining whether it 
is a separable state or not. Of course, for larger dimen- 
sions the subproblem of minimising the respective target 
function becomes increasingly more time consuming. Not 
only the EoF itself, but also an optimal ensemble realis- 
ing the state can be calculated. 

We have extended results on the concurrence and EoF 
of 2 X 2 systems by also including the cardinality of the 
optimal ensembles. More importantly, we have tried to 
generalise the concept of concurrence to general systems, 
and have shown that this generalised concurrence has 
potential to supply a fast test for separability of general 
bipartite states. 

In the future, we will use the presented methods to 
generate more numerical results about EoF of higher- 
dimensional states, for example, to chart the "unknown 
territory" of bound-entangled states, or just as a means 
for testing various conjectures. Furthermore, the varia- 
tional characterisation of EoF could be useful in proving 
or disproving that EoF is additive. Another interesting 
topic for future work is trying to prove the conjectured 
sufhciency of the generalised concurrence test for separa- 
bility. 
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APPENDIX A: PRESELECTION OF T MATRIX 

The topic of this appendix is a method for reducing 
the set of T matrices over which the minimum (^) has to 
be taken in a separabiUty test. In some cases the method 
already yields the optimal T matrix without need for 
performing a minimization procedure. This method is 
based on a method used in blind identification for array 
processing 

Consider the expression 



V B 



p,q 



where Apg'' is as defined in and Bi jis a symmetric 
matrix. When we substitute equation (0) in it, we get, 
using bilinearity of C: 

BpqApq 

p,q 
R 

= ^ BpqC{^/rn^4^ .y/m^cj)'^) 

p,q=l 

R K 

p.q—1 k,l — l 

K 
k,l=l 

Note that, just like B, (T^BT*) is also symmetric. 

Suppose that the state p is indeed a separable one; 
then there exist matrices T leading to a product state 
decomposition, i.e. to C{ip'',ip'') being identically zero for 
k = I. Consider one such T. There exist symmetric 
matrices B for which (T^BT*) is diagonal, say equal to 
some A. Indeed, by right-unitarity of T one just has to 
take 

B = TAT^. (Al) 
Using such a B in the above expression, we find 

Y^Bpq^-^^O, (A2) 



for all a. 

We can now reverse the reasoning and say that any T 
leading to a product state decompositio n ca n be found 
from some symmetric B that satisfies (A_2). That is, 
instead of searching for a T in the complete set of unitary 
mat rices, we o nly have to consider T that follow, using 
( |A1| ) and (A2), from such B. If T is square (that is, 
K = R),T is unitary, and since B = TAT^, 



BB* 



T\A\'^Tl 



Hence, the column vectors of T must be the eigenvectors 
of BB*. Given then all the symmetric matrices B that 
satisfy equation (A2), we only have to consider T ma- 
trices whose column vectors are the eigenvectors of one 
such B. 

We will now show that under some conditions the re- 
duced search space contains nothing but the optimal T, 
so that no search has to be done at all. In that case, one 



just has to take one B satisfying equation (A2), and con- 
struct a T from its eigenvectors. The first requirement 
for this is that the cardinality K must equal the rank R, 
so t hat T is then unitary; the reason is that otherwise 
(Al) has no unique solution. Let us suppose that the 
first P {P < K) statevectors in the ensemble realising 
p are product vectors: \'tp'') = \a'') (g) < k < P. 

Therefore, C{ip'' will be zero for k < P. Now, the 
matrices C{ijj'' ,tjj'') for k < I and k = 1 > P are in gen- 
eral (that is: for all states except for a subset of measure 
zero) linearly independent, as long as the number of ma- 
trices does not exceed the number of matrix elements. 
If the latter requirement is not fulfilled, then of course 
a dependence must exist between the matrices. If the 
requirement is fulfilled then the matrices can still be de- 
pendent provided the K vectors (being m — KNiN^ 
complex variables) satisfy a system oi Ni{Ni — l)N2{N2 — 
\) / A - K{K - I) /2 - K + P + \ polynomial equations 
of degree d — K{K — 1) + 2{K — P) (each equation cor- 
responds to a minor of rank K{K — l)/2 + K — Poia. 
matrix containing {'^'^ S'^°'^'^)ki as elements). Using the 
Schwarz-Zippel theorem we find that the set of vec- 
tors obeying just one of those polynomial equations has 
measure zero with respect to the set of all possible sets of 
K vectors. A fortiori, this also holds for the set of vectors 
obeying all polynomial equations. We thus get a second 
requirement for the automatic optimality of T, namely 
that the cardinality K must satisfy the inequality 

K{K - 1) ^ ^. ^^ N.jN, ~ 1) N2iN2 - 1) ^ ^^3^ 



It then follows that q BpqApq can only be zero if 

(T'fBT*)ki = for all k '^ I and k = I > P. In other 
words: (T^ BT*) is necessarily a diagonal matrix for any 
B satisfying (K3), and any T obeying (Al) is optimal. 

We have not investigated whether this technique for 
reducing the search space is also applicable for calculating 
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the EoF; that is, whether some T that is optimal w.r.t. 
(0) can be found in the reduced search space. 



APPENDIX B: CALCULATION OF THE 
GRADIENT OF THE AVERAGE 
ENTANGLEMENT 

The geodesic on the unitary manifold along a direction 
X (skew-Hermitian matrix) and passing through Tq is 
given by = roexp(eX), or To(l + eX), for small e. 
The gradient of a scalar function on the manifold can be 
calculated from the variation of the function along the 
geodesic, using 



de 



(V/,X) 



To avoid notational clutter, we have set Tq equal to 1 in 
the rest of the appendix. 

Let us recollect that the function of T which is 
to be minimised is g{T) — J2k^i^k(T)), where 
G{A) = -Tr(Alog2(A/Tr(A))) and Afc(T) = 

Lemma 3 For Hermitian A and B, 



G{B,A), 



'=0 



where 



g{B, A)^~ Tt{B log2 A) + Tr(S) log^ Ti(A). 

Proof. We use the following formula from (formula 
6.6.31), which applies for a Hermitian matrix A(t) func- 
tion of a parameter t with eigendecomposition A{t) = 
U{t)A{t)U{ty , and for analytic functions /: 

|/(A(t)) ^ U [(A/(A„ A,))., o U^A'U] C/t. 

Here, o is the Hadamard product and Af(\i(t), Xj{t)) 
are the "divided differences" 



A/(A.(i),A,(i)) 



/(A.(f))-/(A,(f)) r ■ , ■ 
f'iMt)), iOTZ=J. 



For A{t) = A + tB, it follows that: 



^Tr/(AW) 



t=o 



Setting f{x) — h{x) = 
we have f'{x) — —(1 



= ^A/(A,(t),A,(t))(f/ti3C/),, 

i 

= Trf'{K)U^BU 
= Tr/'(A)B. 

-xlog2(a;), so that f{A) ^ H{A), 
In x) / In 2 and 



dt 



TTH{A + tB) 



Tr(l + lnA)B/ln2. 



Furthermore, 
d 



dt 



h{A + tB] 



-(l-|-lnTrA)TrB/ln2, 



so that the lemma follows. □ 
Proceeding in a similar fashion, we can expand Afc(Te) 
up to first order in e. Putting QP"^ — y/fnprr^ij)^ : 



p.q 

p,q 

p 

where we have used the fact that X is skew-Hermitian. 
Inserting this expression in ^ G{Ak{Tf_))\^_^ we see 
that Q'='= serves the role of "A" and J2pC^pkQ'''' - 
XkpQ^P) that of Exploiting linearity of Q with re- 

spect to its first argument, we arrive at the expression 

= XpkigiQp^Q'"') ~ g{Qp\Qn) 



p,k 



(in the last term we have interchanged the indices k and 
p). Therefore, 



^pk /^kk\ 
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FIG. 1. Entanglement of formation for Horodecki states in 
function of a and e; linear scale. 

FIG. 2. Entanglement of formation for Horodecki states in 
function of a and e; logarithmic scale. 

FIG. 3. Entanglement of formation for Horodecki state 
a = 0.225 in function of e; linear scale. 

FIG. 4. Effect of cardinality on calculation of entanglement 
of formation. 



11 



. 0135 
0.013 

. 0125 
0.012 

. 0115 
0.011 



